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Abstract. The Lagrangian theory of structure forma- 
tion in cosmological fluids, restricted to the matter model 
"dust" , provides successful models of large-scale structure 
in the Universe in the laminar regime, i.e., vifhere the fluid 
flow is single-streamed and "dust" -shells are smooth. Be- 
yond the epoch of shell-crossing a qualitatively different 
behavior is expected, since in general anisotropic stresses 
powered by multi-stream forces arise in coUisionless mat- 
ter. In this paper we provide the basic framework for the 
modeling of pressure-supported fluids, restricting atten- 
tion to isotropic stresses and to the cases where pressure 
can be given as a function of the density. We derive the 
governing set of Lagrangian evolution equations and study 
the resulting system using Lagrangian perturbation the- 
ory. We discuss the first-order equations and compare 
them to the Eulerian theory of gravitational instability, 
as well as to the case of plane-symmetric collapse. We ob- 
tain a construction rule that allows to derive first-order 
solutions of the Lagrangian theory from known first-order 
solutions of the Eulerian theory and so extend Zel'dovich's 
extrapolation idea into the multi-streamed regime. These 
solutions can be used to generalize current structure for- 
mation models in the spirit of the "adhesion approxima- 
tion" . 

Key words: Gravitation; Hydrodynamics; Instabilities; 
Methods: analytical; Cosmology: theory; large-scale struc- 
ture of Universe 



1. Introduction 

Zel'dovich (1970, 1973) initiated the use of Lagrangian 
coordinates for the construction of analytical models of 
cosmic structure formation. He pointed out that the so- 
lution of a force-free continuum of "dust" (i.e. pressure- 
less and non-gravitating matter), given in terms of tra- 
jectories of continuum elements (see, e.g., Zel'dovich & 
Myshkis 1973), can be rescaled to give the correct solu- 
tion for gravitational instability in the linear limit. He 
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suggested to follow those trajectories further up to the 
epoch when shell-crossing singularities in the fluid flow 
develop. Doroshkevich et al. (1973) confirmed the valid- 
ity of Zel'dovich's approximation also within the full self- 
gravitating system of equations. To be more precise they 
performed a self-consistency check in one of the four New- 
tonian field equations (see Buchert 1989 for a discussion 
and extension of this self-consistency argument). It was 
also known (Zentsova & Chernin 1980) that Zel'dovich's 
ansatz provides a subclass of exact solutions in the case of 
plane-symmetric collapse on a Friedmann background. In 
(Buchert 1989) it was then shown (going back to an ear- 
lier work on the fully Lagrangian formulation of the basic 
equations by Buchert & Gotz 1987) that his ansatz also 
provides a subclass of exact 3D solutions for a special class 
of initial gravitational potentials that are composed of sur- 
faces of vanishing Gaussian curvature. These solutions fea- 
ture maximally anisotropic, locally one-dimensional col- 
lapse supporting Zel'dovich's original discussion of the for- 
mation of "pancakes" as collections of volume elements 
which degenerate into small sheets and later enclose a 
three-stream system of the flow. These findings may be 
restated in the framework of a Lagrangian perturbation 
theory in which Zel'dovich's model is recovered as a sub- 
class of the first-order solutions (Buchert 1992; for reviews 
see Bouchet et al. 1995, Buchert 1996, Sahni & Coles 1996, 
Ehlers & Buchert 1997 and references therein). 

The formation of singularities (Arnol'd et al. 1982) is a 
consequence and drawback of the matter model which al- 
lows fiuid elements to overtake freely. On the other hand, 
N-body simulations show (e.g. Doroshkevich et al. 1980) 
that N-stream systems create additional multi-stream 
forces which hinder the central fraction of fluid elements to 
escape from high-density regions; as a result they fall back 
onto the central region performing oscillations about the 
center: an internal hierarchy of nested caustics is formed 
that continues to create new structures down to smaller 
scales (Gurevich & Zybin 1995). 

A phenomenological model based on Burgers' equation 
was then suggested to overcome the drawback of "dust" 
models by inventing a Laplacian forcing that acts in the 
desired way (Gurbatov et al. 1989). This model features 
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"adhesion" of already formed structures and so models the 
action of self-gravity of multi-stream systems. A recent 
suggestion of how to derive this "adhesion approximation" 
from kinetic theory of self-gravitating coUisionless sys- 
tems shows that it is indeed possible to obtain adhesion- 
type models by taking into account multi-stream forces 
(Buchert & Dommguez 1998). However, these models form 
a wider class in the sense that multi-stream forces are, in 
general, anisotropic unlike in the case of Burgers' equa- 
tion. Moreover, equations involving a Laplacian are ex- 
pected for the gravitational field strength and not, as is 
the case for the "adhesion model" , for the velocity. 

Inspitc of the insight gained by deriving the "adhesion 
model" from first principles, we shall, in the present work, 
adopt the restricting assumptions that the multi-stream 
force acts isotropically (in the sense of an idealization), 
and that the dynamical pressure due to multi-stream sys- 
tems can be represented as a function of the density. For 
these assumptions we derive the general Lagrangian evo- 
lution equations and expand them to first order using the 
Lagrangian perturbation approach. 

A comment concerning the isotropic approximation of 
the multi- stream force is in order: our investigation covers 
pressure-supported fluids as they are studied in hydrody- 
namics. Here we do not restrict ourselves to the descrip- 
tion of matter models which are supported by the usual 
thermodynamical pressure (i.e., by short-ranged interac- 
tions). Rather we follow an approach that we consider to 
be relevant for coUisionless systems: a dynamical pressure 
force arises due to the action of a "multi-dust" region on 
the bulk motion. As remarked above such a forcing is in 
general anisotropic and the evolution equations govern- 
ing the anisotropic stresses involve all velocity moments 
of the distribution function (Buchert fc Dommguez 1998). 
Hence, a closure condition is needed which, especially for 
self-gravitating fluids (e.g. Bertschinger 1993), is not obvi- 
ous and cannot be formulated along the lines of collisional 
systems. At the current stage of knowledge about this sub- 
tle problem we advance the phenomenological closure con- 
dition that multi-stream stresses may be idealized by a 
scalar function p that is given in terms of the density of 
the continuum. This opens the possibility to immediately 
access the multi-streamed regime phenomenologically. For 
example, Buchert & Dommguez (1998) found that the as- 
sumption of small velocity dispersion singles out a rela- 
tionship p oc Q^^^ that is solely based on the dynami- 
cal equations; the matter model p (x g'^ yields, together 
with the restriction to parallelity of peculiar-velocity and 
-acceleration, the "adhesion approximation" (Gurbatov et 
al. 1989). Other relationships are put into perspective by 
Buchert et al. (1999); particular matter models may lead 
to soliton states which, due to their persistence in time, 
could dominate the architecture of cosmic structure. 

Although the isotropic idealization is simple-minded, it 
may actually be a good approximation in practice, as has 



been recently found for high-density regions in numerical 
simulations (Colombi, priv.comm.). 

We proceed as follows: In Sect. 2 we develop the La- 
grangian theory for pressure-supported fluids and illumi- 
nate the structure of the basic equations by the plane- 
symmetric case; in Sect. 3 we then move to some detailed 
studies of the basic equations: a perturbative treatment al- 
lows to derive a first-order evolution equation for fluid dis- 
placements. Note that this linear equation in Lagrangian 
space embodies Eulerian nonlinearities. This allows us to 
formulate an extrapolation rule that extends solutions of 
the Eulerian theory of gravitational instabilities into the 
nonlinear and multi-streamed regime. This rule rests on 
the spirit of Zel'dovich's extrapolation idea that led to the 
celebrated "Zel'dovich approximation" (Zel'dovich 1970, 
1973). We also discuss the special case where the gravita- 
tional collapse is plane-symmetric on an isotropic back- 
ground cosmology. We summarize the results and provide 
an outlook in Sect. 4. 

2. The Lagrange— Newton System with pressure 

2.1. The Euler-Newton-System with pressure 

Eulerian coordinates are the most widely employed ones 
in physics. Here we use x to define the position that a 
fluid element occupies at the moment < in a non-rotating 
Eulerian coordinate system. If we consider x and t as inde- 
pendent variables, the total time derivative of any vector 
field a = a{x,t) equals: 

da da , 

- = - + (.. V)a . (1) 

The field v{x,t) corresponds to the mean velocity in a 
kinetic picture. Eulerian coordinates are used to describe 
fields, e.g. the spatial mass distribution g{x,t). In this 
description we are not interested in the fiuid elements that 
produce this density and we do not know their trajectories. 
In Subsect.2.2 we introduce Lagrangian coordinates which 
convey this information. 

The basic system of equations describing a self- 
gravitating medium is the so-called Euler-Ncwton- 
System. We add the accelerative term — which takes 
into account the isotropic part of the velocity dispersion 
tensor. We have the following evolution equations for the 
mass density g and the velocity v {b denotes the acceler- 
ation field): 

^ = -V-igv) , (2a) 

dv 

= —(v ■^)v + b , where b := g . (2b) 

at Q 

The gravitational field strength g is a solution of the (New- 
tonian) field equations: 

V X g = , (2c) 

V • g = A - AirGg , (2d) 
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where A is the cosmological constant and G the Newto- 
nian gravitational constant. The evolution equation for g 
represents the conservation of mass (continuity equation). 
The second equation (for v) is Euler's equation and de- 
scribes the balance of momentum. In a kinetic picture (2b) 
is known as Jeans' equation. 

To complete the Euler-Newton-System with pressure 
we need an "equation of state" that relates p with the 
dynamical variables g and v: we employ the assumption 
p = a{g) with a' := Hence, the field equations become 



Vp 



= Vxg = Vxb + Vx 

Q 



= V X b + -V X [Vp) + ^Vp X Vg 
g 

= V X b : 



(3a) 



A - 47rG£> = V ■g = V -b + V ■ — 

Q 

= V • 6 + Vp • V- + - Ap 

g g 

^V-b+ia"-^)^^ + ^Ag . (3b) 
g g g 

The equation for the curl of g does not change with the 
additional accelerative term for that particular class of 
"equations of state" . 



2.2. Transformation tools 

Lagrangian coordinates are assigned to fluid elements and 
do not change along the flow lines: we choose the positions 
X of the elements at a time to as coordinates {X — x\tg); 
thus we have indexed each fluid element. If we use X and 
t as independent variables, the total time derivative be- 
comes ^ = ^1^' i-^-' convective term disappears in 
this description. In order to connect the Lagrangian coor- 
dinates with the Eulerian ones, we introduce the position 
vector field /, 



x^f{X,t), where X := f{X,to) 



(4) 



For later convenience we introduce the symbols " , " 
and " I " in order to distinguish Eulerian and Lagrangian 
differentiation, e.g. =: Va.t and =: fa\b- V denotes 
the Eulerian and Vo the Lagrangian nabla operator. In 
the following we also apply the summation convention. 

The Jacobian matrix, which reflects the deformation 
of a volume element, is of key importance for the trans- 
formation of coordinates. Here we use the notation: 



J := det Jij 



det/,| 



(5) 



In order to change the Eulerian differentiation to one with 
respect to Lagrangian coordinates, we have to deflne the 
inverse transformation of (^): 

X = h{x,t) , h = f~^ . 



Remembering the definition of the adjoint matrix, 
adJjTc = ■^£kim£jpqfp\ifq\m, we get 

^fcj — Jjk ~ 'J^'iJjk — -^£klm£jpqfp\lfq\m 'j (6) 

Skirn denotes the totally antisymmetric Levi-Civita tensor; 
note the useful relation £ijk£iim = &]i5km - SjmSki {5ij = 
1 for z = j and 5ij = else). In order to shorten the 
equations it is convenient to use functional determinants: 



J{A,B,C) 



d{A,B,C) 



£klmA\kB\iC\ 

m ■ 



dA 


dB 


ac 


dXi 


dXi 


dXi 


dA 


dB 


dC 


0X2 


8X2 


0X2 


dA 


dB 


dC 


dXs 


0X3 


0X3 



(7) 



Because of the properties of the Levi-Civita tensor, or 
the definition of a determinant as a multilinear and alter- 
nating map, respectively, we have some tools that make 
work easier, e.g.: 

J{A + D,B,C) = JiA,B,C)+J{D,B,C) 
J{A,B,C) = -J{A,C,B) 
J{A,A,C) = 
J{A ■D,B,C) ^ D- J{A, B,C)+A- J{D, B, C) . 

Therefore, we can transform any tensor j as follows: 

O^iJ — ^i\k^k,j — ^i\kl2j£klm£jpqfp\lfq\m 



(8) 



For further details on Lagrangian evolution equations the 
reader may consult the introductory tutorial of Buchert 
(1996) and the review by Ehlers & Buchert (1997). 

2.3. The Lagrange-Newton-System with pressure 

Starting with the Euler-Newton-System (with pressure) 
we obtain the Lagrange-Newton-System with the help of 
the transformation x = f{X, t). 

As the continuity equation ( pa| ) represents the conser- 
vation of mass M — Jj^^ gclc^x ~ gJclc^X within any 
domain Dt that may change in time, i.e.. 



dM 
"dT 



O^A/' gJd'X^f ^(gJ)d'X 



Dt 



J 



d-'x, 



we conclude with go := g{X,t()), Jq :~ J{X,t()): 

gj =: C{X) = goJo . (9) 



For our definition of Lagrangian coordinates (|4|) we have 

dXi 



Jq — det ■$^'1 to — 1; nevertheless we will use C{X) for 
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the time being in order to allow for a later relabelling of This system of equations can be closed with the help of 

trajectories. the exact integral (^. 

According to the definition of / the Eulerian equation 

(pbl) reads: ^ , , ^ ^. 

' — ' Z.4. Annotations 

b = f{X, t) . (10) g^^jg^ii-^^-j^g ^jjg integral (|) into (p^) we obtain: 



As we have the exact integrals (|9|, |10|), g and v are 
no longer dynamical variables in the Lagrangian picture; ^ _ ^^q^ _ _}_£ . f f ) 

/ attains the status of the only dynamical variable, if we 

demand that p is a functional of /, which is obviously true , o;' ^/ 't//^ f f \ f f \ 

for the assumption p = a(£i). 2CJ"^ 

Applying the first integral (^) and the transformation 
of the Eulerian differentiation m) to the equations for 2 ' 



o 73 '-^ ("-^ ■^P' ■^9)' •^P' •^9) 



g (|2q, 2d) we obtain the Lagrange-Newton-System with 1 f n J 



pressure, i.e. a set of four coupled nonlinear partial differ- ^ 2CJ3^^^' ^'^^'^^^^ -^P' -^^-^ (^"^ (^"^ 

ential equations for the trajectory field /. We first obtain: / j 

- ^^(^' /P' h)^iJ^ /P' /«) (2a" + ^a' 



= £/,.jJ»j + -nrQ,^P,l (13) 

= -jJifi, ft, fh) Let us now introduce the following approximation which 

^ is common in cosmology: consider the definition g{X) ~ 

+ 2j2 2 ^3P'?^(^' •^'i' /p' /«) ' (11^) £'h(1 + 5(^)), where gin denotes the background density of 

^ the mean matter distribution and S the density contrast. 

A - AnGg = V ■ g = V ■ / + V • — At the time i = to we have ga{X) = gHo(l + Sa{X)). 

^ As the density contrast is numerically very small at the 

_ . ^^(j f f ) time of recombination (i.e. to = tree), go — gw is a useful 

2>/ approximation. Therefore we get C\i — {goJo)\i —■ Cii\i — 

_l_ J[ ^ c kiJjp Ik fi) f f ) ^ ^^'^^ alternative exact argumentation see Appendix 

we use Ch instead of C furtheron to emphasize 

_ 1 Ji f f f ) homogeneity. The Lagrange-Newton-System then may 

- ^j^^Pi^'-J^'-'P'-'^' be written in a much simpler form, since all functional 



1 ^/ £ £ \ rr/ T J- J- \ determinants with C disappear: 

J[P,Jp,Jq)J[J-JpJq) 



2J3g 



1 r r.^f f JifrJrJh) = , (14a) 



A-47rG^ - —ejpqJif,Jp,.fq) 



1 I 2 J 

+ 7n2-.JiJ (p^fp^f,) Jp J,)- (lib) 



In the special case p — a{g) we get our final set of equa- ^ / 2 J 

tions: + i^J{J, fp, fq)J{J, fp, fq) + 

= (V X /);, = ShJ^k, = jj{f\, f,,h) , (12a) (14b) 

A — AwGg = ——£jpqJ{fj,fp,fq) Looking at Eq. ( |l2b| ) we conclude that a (not necess arily 

physically relevant) possibility that simplifies Eq. ( |l4b| ) 
-J{g, fp, fq)J{J, fp, fq) further is the following: 



2J 
a' 

2J^g' 
a' 

a 



2J 2 1 

J{Q,fp,fq)J{Q,fp,fq) a" ^ -—a' = — a' , i.e. , p = a- + C2 .(15) 



2j2 ^ (fi*' -A" /?)' /p' /?) This "equation of state" plays a special role as will become 

„ clear in Subsect.3.2. 

" -Ji0,fp,fq)Ji0,fp,fq) ■ (12b) 



2J^g 
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2.5. Reduction to planar symmetry 



In order to learn more about the structure of Eqs. (12a 



12b ) it is useful to look at restrictions in symmetry: we 
restrict x — f{X,t) to planar symmetry, i.e.,[] 



(16) 



the only direction in which motion takes place is the Xi- 
direction; /2 and /a are constant. Consequently, there 
are no velocities and accelerations in the X2- and X3- 
directions (in particular (72 = 53 — 0). So, J simplifies 
to 

dxi 



J = 



dXi 



The four field equations that we have in 3D are reduced 
to the equation for the divergence of gi , 

A - 47rGf? = A - 47rG-^ = g,., = ^ ; 
/i|i 

the curl of g vanishes identically. With the help of 

gi = oi + — = /i + —r- = Ji + ■ 



we get 



f/iii 



A/i|i -47rGeo = All + (^)|i 



Qo 



Let us define Gi|i := A - AttGqo {G{Xi) := g{Xi,to)) in 
order to cast the last equation into a total divergence, 



(/i-A/i-Gi + AXi + ^)|i = 

90 



hence. 



/i - A/i = Gi - AXi 



Ell 

eo 



neglecting an irrelevant function of time. With the relation 
P = ctis) the pressure-force term becomes 



P\i a' 



Qo 



Qo 



Qo V/i| 



Giiii 



/i|ii 

(Gill - A) An "727 



Altogether we arrive at 



/i - A/i = Gi - AXi + a- 



Gi 



111 



/2 (Gi|i-A)/i 



(17) 



or (with ignorance of dimensions) 

A - A/i = Gi - AXi + fin f 77^^^ 

V \Gi|i - A 

(18) 

For the case A = and a' = const. Gotz (1988) has shown 
that (^8|) can be mapped to the Sine-Gordon equation, 
which is a well-studied equation that admits soliton solu- 
tions. 



^ In this subsection we use Jo = 1, hence C(X) = Qo{X) 



3. Lagrangian perturbation approach 

3.1. The perturbation ansatz and linearization 

In standard cosmology we invoke a homogeneous defor- 
mation of the continuum: /hi = aij{t)Xj that is isotropic: 
aij{t) = a{t)6ij; thus, with /h = a{t)X and the homoge- 
neous density Qh '■= ^ applied to the Lagrange-Newton- 
System with pressure we are left with 



A-47rG 



Gh 



(19) 



a 

The first integral of this equation is Friedmann's equation: 
^2 -L const. SttG^ + A 



a 



(20) 



In order to describe structure formation we consider small 
{0{e)) deviations p{X,t) from this homogeneous and 
isotropic motion; we use the ansatz 



f{X,t) - a{t)X+p{X,t) 



(21) 



and suppose that perturbation theory is justifiable. The 
assumption that the perturbations are smooth leads us to 
\Pi\j\ ^ a and \pi\j \ ^ d. Inserting our ansatz (21) into 
the first equations (14a) we have 



= (Vxg)„ = J(/„/„/,0 

= J{aXi + Pi, aXi + Pi, aXh + Ph) 

= daJiXi,pi,Xh) + aJ{Xi,pi,ph) + a^Jipi,Xi,Xh) 

+ aJipi,X^,ph) + aJ{pi,pi,Xh) + Jipi,Pi,Ph) ■ (22) 

In the first-order approximation we neglect quadratic and 
cubic terms, hence: 

daJ{Xi,pi,Xh) + a'^J{pt,Xi,Xh) = or equivalently, 



Vo X p - -p = 



(23a) 



Linearizing the equation for the divergence of g we 
start with (|l4b| ), i.e., with the assumption G = Gh (com- 
pare Appendix A) . The following terms appear during the 
calculation: 



J 
1 
J 



1 

£jpq >J{fjTfp:fq) 
{Ji fpT fq) 
3 (iT" (</j /p, /g)j /g) 

J{JJpfq)J{JJpJq) 



+a'^p,l, + 0{e^) 
1 

a3(l + ip,|,) + 0(e2) 



1 



1 



3-(l--p,|,) + 0(£7 
a-'" a 

6aa^ -I- Aaapiii + 2a^Pi\., + 0{e^) 

,4„ I /0/',2\ 



£apqa Pt\ia + 0{e 
2aV.K„- + O(e') 
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(remember: SabcSabc = 6 and EabcSabd = 2(5cd)- In the end 
we find: 

„Ch „C'h a a 1 
= -A + AttG^ ~ 47rG— 3-p,|i + 3 + -p,\i 



(23b) 



(The calculation without using C = Ch may be found in 
(Adler 1998).) In comparison to the matter model "dust" 
(Buchert 1992) there is only one additional term ^Vi\ijj — 
^Ao(pi|i). Inserting the zero-order equation ( |l9|) (thus 
neglecting backreaction of the first-order displacements 
on the background solution) yields 

^ , „Ch a 1 .. a' 

= -47rG— + -pi|i j,Vi\%n ■ 

(JL G/ Gj Gj 

Introducing the comoving displacement field P :— 
representing the inhomogeneous deformation that is scaled 
with the expansion, we get: 

= -^-KG^Piu + Piu + 2-Pju j-^lw ■ 

This equation can be written as a total divergence: 



Vo • ( P + 2-P - 47rG^P - ^AoP ) = 



Writing Equation (23a) also in terms of the scaled dis- 



placement field, and also inserting the zero-order solution, 
we arrive at our main result: 

For vanishing harmonic parts (see Ehlers & Buchert 1997) 
we obtain the following final set of equations that com- 
prises the linear perturbation theory for the scaled dis- 
placement field P = P^ + P^ ( Vo X P^ = ; Vo • P^ = 
0): 



PT ^1-PT 

a 

pL + 2-Pi-47rG^P^ 



^AoP- 



(24a) 
(24b) 



Remark: 

Remember that a' ~ Ci {o) where g depends on the defor- 
mation P; thus, the r.h.s. of Eq. (24b) is not linear in P^. 
So, for given p = ol{q), a'AoP^ has to be linearized with 
respect to P^ in the context of linear perturbation theory, 
i.e. a' = const. This fact amounts to a major difference 
compared to the "dust" case: in the latter, the dependence 
on Q completely disappears in the equations and the ex- 
act density may be obtained for any perturbative solution. 
Here, the density is implicitly linearized and this fact has 
already been used during the above derivation. 



3.2. Comparison with plane-symmetric inhomogeneities 
on a Friedmann background universe 

Let us start with trajectories, which are scaled with the 
expansion F = Reduction to planar symmetry now 
means: 

Fi=Xi+ Pi{Xi,t) ,F2 = X2,F3^ X3 . (25) 
The functional determinant J reduces to 

J = a3(l + Pi|i) , 
and, thus, J is independent of the coordinates X2 und X^. 
With the simplification C = Gh we derive from ( |l4b| ) 

Gh 2a^'d, 



(a(H-Pi|i) + 2aPi|i +aA|i) 



I.e., 

A - 47rG 
„ d 



J3 

Gh 



Gh 
7^' 



2J , 



Gh 



a )a J|i J|i 



fl3(l + Pi|i) 



a 



All 



All 



a 



a a(l + Pi|i) 1 + Piii a^l + P^,r 
Ch . 2(1 + Pi|i)a3 , 



A I 



a-^(l + Pi|i)' 



Gh 



-a')Pi|iiPi|ii- 



Using the zero-order equation (19) and multiplying with 
(1-KPiii) yields 



Pi + 2-Pi - 47rG^Pi 



a 

Gh 



,1 a'(l + Pi|i)- 



^A|i 



11 



2(1 + Pi|i)a3 , 



a')Pi|iiPi|ii 



a^(l + A|i)^' Ch 
The r.h.s. of this equation can be written as a total deriva- 
tive; hence we get for vanishing harmonic parts. 



Pi + 2-Pi - 47rG^Pi 



a2(l + Pi|i) 



2P1I11 



(26) 



By linearizing this last equation we recover Eq. ( 24b ) for 
one-dimensional inhomogeneities. Writing the coefficient 
on the r.h.s. as a function of g, 

a'{g) _ a'{g) ( ^ ^ 



x{q) 



a2(l + Pi|i)2 a2 

we see that we get a linear equation in the case of the 
special "equation of state" (p^. 

Let us finally restrict our problem to the limit of no 
expansion a{t) — a{to). We have gu ~ Gh and /i = 
a{to){Xi+Pi) from (|f|). Therefore, 



fi - 4^G 



Gh 



^3-(/i-a(io)Xi) 
If we additionally appreciate that (|l9|) now reads A 



a[to)^ /i|i 



A-kGCh and Gi|i = A - 47rG£iH0 = 0, we recover Eq. ( |18| ) 
in the limit of the static Einstein universe. 
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3.3. Comparison with Eulerian linear theory 



4. Discussion of Results and Outlook 



It is a standard procedure in cosmology to scale the Eule- 
rian coordinates with the homogeneous and isotropic ex- 
pansion (dh = Hx with Hubble's function H :— fj, i.e. 
to use the following transformation of coordinates 

q = Q{x,t) = ^ . (27) 

The comoving differential operators (Vq)i 
^ \q become 
1 



and 



d_ 

di 



- ^ 
^ dt 



The deviations from the homogeneous fields are called 



u{q,t) 
w{q,t) 

6{q,t) 
where vh 



v{q,t) 
g{q,t) 



- Quit) 



Quit) 
aq and gjj 



peculiar velocity, 
peculiar acceleration and 

density contrast, 

= aq denote the Hubble- 



velocity and the Hubble-acceleration. We obtain the 
Euler-Newton-System with pressure for the peculiar 
quantities by making use of the homogeneous solutions 
(see Peebles 1980): 



-{l+6)Vq-U = 



Hu 



Vq X W 
V„ ■ W 



1 VqP 



(28a) 

' (28b) 

a £ih(1 + o) 

(28c) 

= -AirGgYiSa . (28d) 

In the following we assume that all peculiar quantities 
and their derivatives are small (C'(e)). Using perturbation 
theory to the first order in e we get the linearized form of 
the equations: 



d_ 

dt 



u = 0{e^ 



d_ 

dt 



^ q 
V 



Hu = w 
X w = 



1 VqP 

a gu 



'AnGgnSa 



(29a) 

(29b) 

(29c) 
(29d) 



The linearized equation for the evolution of S can be cal- 
culated by applying — ^Vq to the second equation and 
inserting the first and forth equations: 



dt^ 



d 

S + 2R- 
J at 



S - AnGgaS = 



a^gu 



AqS 



(30) 



/ dp 



is defined 



where we have set p = a{g), and a' = 
in terms of the "speed of sound" Cs as is usual in a hydro- 
dynamical medium. Here, we have a' = const., i.e., p oc g, 
in order to obtain a linear equation in S. 



In the last section we have derived three different equa- 
tions: 

— From the linearization of the Lagrange-Newton- 
System with pressure we have for the longitudinal part 
of the inhomogeneous deformation P{X,t) = ^^^j^: 



pL 



AnGguP^ 



(31A) 



remember gn '■= 

The plane-symmetric solutions on a Friedmannian 
background obey: 



Pi + 2-Pi - AnGgnPi = 
a 



a + Pill)'- 



:Pl\ 



(31B) 

The Eulerian linear theory in comoving coordinates 
q — leads to a partial differential equation for the 
density contrast d{q,t): 



a d 
q a ot 



S - AirGgnS = 



(31C) 



Acco rding to the equivalence of the equations (|31A|) and 



( t^lCt ) up to the time derivative operators ^ — -Six and 



■§[\q, we see that, with the help of the already known re- 
sults for the density contrast in the Eulerian linear theory 
(e.g. Haubold et al. 1991), solutions for the inhomoge- 
neous deformation P^ in the Lagrangian linear theory 
can be constructed. In the case p — a. class of exact 
3D solutions has been found with the help of this method 
(see: Buchert 1989). But, in contrast to the case p — 0, 
this is not to be expected here, since extrapolation of the 
solution of the linearized equations does not yield an ex- 
act solution of the planar problem as well; this is easy to 
see from Eq. (31B): the pressure term produces (except 
in the special case (p^)) nonlinear terms in P already for 
plane-symmetric inhomogeneities. 

We emphasize the special role played by the "equation 
of state" p = ^ + 02. 

For further applications we can proceed as follows: 
We look at solutions S^{q,t) of the differential equation 



(31C) and use the instruction q 



X, 5' 



Pt to 



const ruct solutions Pl'{X,t) of the differential equations 
(31A). As the Lagrangian description implicitly respects 
nonlinearities, our construction rule allows to build non- 
linear models of structure formation that take dynami- 
cal pressure forces into account. Since solutions to Equa- 
tion ( ^1A| ) constitute the Lagrangian extrapolation of the 
Eulerian linear theory, they are built in the same spirit 
as Zel'dovich's approximation in the case p = Q (see: 
Zel'dovich 1970, 1973; Buchert 1989). Viewed together 
with the derivation of the "adhesion approximation" by 
Buchert & Dommguez (1998), these solutions may be used 
as first approximations to adhesive gravitational cluster- 
ing in the weakly nonlinear regime. 
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Appendix A: Alternative argument for the approx- 
imation C\i — 

We want to give an alternative argumentation for the ap- 
proximation C\i = {goJo)\i = 0- Let us introduce new 
curvilinear coordinates Y by 



Y := A{X) = X + *(X) 



(32) 



Below we shall consider (initially) small deviations '4' from 
X, restricting our argument to the linear approximation. 
We choose the new coordinates in such a way that 

CiY) = go{Y)JoiY) = const. =:Cu ; = Qm ■ 

Thus, taking the conservation of mass into account, we 
have 

dmo = go{X)d^X = gd^Y . 
Writing go{X) = ^)Ho(l + ^o(-^)), the map A is defined 

by 

(l + 5oW)det(^|^) =1 . 
Then we get for — 0{e) and small derivatives: 



l + .o(X).do,(||).dc.(.., + 2|t 



= l + Vo-*(X) + C)(e2) . 

Thus, to first order in e we have: 

5o(X) = Vo • *(X) , 

the particles are initially displaced according to the den- 
sity contrast. With Vo • W — -47rG£<HO(5oa(to) (( ^9d| ) at 
the time t = to) is defined with the help of the peculiar- 
acceleration w at time to : 



1 



4:TTGguoa{to) 



Vo-W, W:=w{X,ta) 



(33) 



To enforce the property = 0, we simply have to relabel 
the trajectories: 

X — > y, x = f{Y, t) with x\t,, = Y = X + *(X) . 

Note that Zel'dovich's original discussion of his approxi- 
mation (Zel'dovich 1970, 1973) as well as subsequent work 
were all using the coordinates Y . We stress, however, that 
the above argument also involves an approximation and is 
only consistent within the first-order solutions; compare 
the discussion in (Buchert 1989). 
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Appendix B: Transformed vector— identities 

The following equations comprise a useful coUeetion of 
formulas; similar expressions may arise in calculations em- 
ploying Lagrangian coordinates. We list these transformed 
vector-identities here, because they may be helpful for fur- 
ther considerations. 

= {V x{Vk))h = jJ{ki,fiJh) 

— ^ipg'-^i'-^i^^ fp,fq)i fii fh) 

~ 2j3 ^ipq'^ (^) fpj fq)^ ("^j fii fh) ! 

= V-(Vxr) = (ij(Tp,/p, /,)),, 

~ 2j2 ^ihqJ^ {J^ (-^P) /p) /g)) fi> fh) 

~ '^j^^ihq>J{Tp, fp, fq)J^{J, fi, fh) ', 



( jShji<J (^r; fri fi)),j 



(r, T^OPq'^i'^ji fpi fq)),h 



+ (TTl^iPI'^^'^h, fp, fq))A 



= (V X (V X T) - V(V • T) + AT)^ 
1 

1 

'2J' 
1 

~ 'j2'^{'^{'^ri frj fi)i fit fh) 

- -j^J{Tr,fr,fi)J{J,fi,fh) 

~ '^j2^hrt£jpqJ {J {Tj , fp, fq), fr, ft) 
+ ^hrtSjpqS {J, /r; ft)<J {Tj , fp, fg) 

+ -^JiJ{Th, fr, ft), fr, ft) 

- ^J{J,fr,ft)J{n,fr,ft) ■ 



